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Abstract 

We present a method to combine next-to-leading order (NLO) matrix ele- 
ments in QCD with leading logarithmic parton showers by applying a suitably 
modified version of the phase-space-slicing method. The method consists of 
subsuming the NLO corrections into a scale-dependent phase-space-slicing pa- 
rameter, which is then automatically adjusted to cut out the leading order, 
virtual, soft and collinear contributions in the matrix element calculation. In 
this way a positive NLO weight is obtained, which can be redistributed by a 
parton shower algortihm. As an example, we display the method for single-jet 
inclusive cross sections at 0{a s ) in electroproduction. We numerically com- 
pare the modified version of the phase-space-slicing method with the standard 
approach and find very good agreement on the percent level. 
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I. INTRODUCTION 



Much progress has been made in the last years in measuring the hadronic final state in 
eP-scattering at HERA with high precision (see [[jj for a recent review) . The theoretical tools 
which are at hand to describe the hadronic final state are basically fixed order perturbative 
calculations, which for most processes are available at next-to-leading order (NLO), or a 
combination of leading order (LO) matrix elements with parton showers (PS), mostly at 
leading logarithmic accuracy, which are implemented in event generators (see [0 for an 
overview of available Monte Carlo programs for both approaches). However, it appears 
that the theoretical calculations are a considerable source of errors in determining physical 
parameters. As an example, the statistical errors in a recent determination of a s from 
dijet production 0] are at the one-percent level. The systematical and theoretical errors, 
on the other hand, are considerably larger and both lie around five percent. Therefore, an 
improvement of the theoretical tools is needed. 

The fixed higher order and the PS approaches have complimentary strengths. The fixed 
higher order calculations reduce uncertainties due to unphysical renormalization and fac- 
torization scale dependences. Wide-angle emission of partons, where interference effects 
between a large number of diagrams may be important, is described well. The PS on the 
other hand allows a description of the cross section in regions where a s becomes large, es- 
pecially in the region of collinear particle emission, by means of a resummation of large 
logarithmic terms. This allows to e.g. describe reasonably well the substructure of jets. In 
addition, the PS can be terminated at some small scale Qo, which allows to attach some 
kind of hadronization model, as, e.g., the Lund model Q, to describe the non-perturbative 
region. It is desirable to find a way of combining the advantages of both approaches into a 
NLO event generator. In a complete NLO event generator one certainly would like to include 
not only the matrix elements at NLO, but likewise the PS in next-to- leading log accuracy. 
Two problems in combining parton showers and matrix elements can be identified. First, 
one has to avoid double counting of events which are included both in the matrix element 
calculations and in the PS. Second, negative weights can occur in the calculation of the 
matrix elements at NLO. Although they are not a principle problem, negative weigths can 
make it in practice difficult to obtain numerically stable results, when they are used as a 
starting point for the PS. 

Basically two strategies can be adopted to combine matrix elements and PS's. Either 
the phase space is split into two parts, using the matrix element cross sections in one region 
and the PS in the other 0-0, or the PS algorithm is modified as to reproduce the matrix 
element cross section in the hard limit ||. In current event generators such as e.g. PYTHIA || 
or HERWIG TO], leading log parton showers are combined with leading order matrix elements 



in the wide angle scattering region. Recently attempts where made to improve the leading 



order accuracy in these approaches by going to next-to-leading order. In |TTJ both, the 
PS and the matrix elements were modified to obtain a smooth merge of the PS to the 
higher order calculation. Collins |T3 suggested a procedure to subtract from the matrix 
elements those parts which are included already in the PS. Since the PS describes the soft 
and collinear region, the divergences associated with the matrix element calculation in this 



region are avoided. Finally, Sjostrand and Friberg [13] suggested an improvement of the 



splitting procedure, by using the NLO matrix elements to calculate the weight for the PS 
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region instead of the LO weight. They suggested to introduce a function which approximates 
the weight in the PS region in such a way that the negative weight problem is largely avoided. 
However, Sjostrand and Friberg only gave an outline of their method without providing the 
approximate function for any specific process, and hence they also did not give numerical 
results. 

In this paper, we take the suggestions in |13[ as a starting point and give further details 
for a general method to get the correct NLO normalization of the PS region. To avoid 
negative weights, we rely on an older suggestion by Baer and Reno |14[, based on the 
phase-space-slicing method. Baer and Reno suggested to introduce a phenomenologically 
determined fixed slicing parameter such that the sum of the Born, virtual, soft and collinear 
contributions is approximately zero. In this paper, we calculate a scale dependent cut-off 
function, which provides a cut-off for each phase-space point, such that the sum is exactly 
zero and the improved scale and scheme dependence of the NLO calculations is preserved. 
The NLO corrections are given by the NLO hard tree-level matrix elements integrated down 
to the cut-off. Since these contributions are positive definite, the NLO weight obtained from 
this method can be directly redistributed according to a PS algorithm. 

The paper consists of two main parts. In the first part we discuss the general method 
to combine matrix elements at NLO with PS's by use of the described cut-off function. In 
the second part we consider as a physically important and interesting example single-jet 
production in deep-inelastic scattering (DIS). Here, the LO contribution is of order a° and 
the NLO corrections are of order a s . Including higher order matrix elements for this process 
will become especially important for diffractive DIS, since the gluon density in the proton 
at small Xbj is large compared to the quark density. Therefore, the photon- gluon fusion 
process, which is a NLO correction to the first order quark-parton model contribution, will 
not necessarily be small. We explicitly construct and numerically study the cut-off function 
for this case. Finally, we summarize our results and give an outlook for future developments. 



II. GENERAL METHOD 

A. Jet cross sections at next-to-leading order 

We start by summarizing the procedure for the numerical evaluation of an inclusive n-jet 
cross section in NLO QCD. The first step is to select a jet algorithm, which defines how 
partons are recombined to give jets. In the following we take for definiteness the invariant 
mass Sij of two partons i and j and define the n-jet region such that Sij < s min , with some 
kind of minimum mass s min and likewise the (n + l)-jet region such that > s min for all i, j. 
The LO process for the production of n jets consists of n final state partons and obviously 
does not depend on the jet definition. This dependence only comes in at NLO. The 0(a s ) 
corrections to this process are given by the ultraviolet (UV) and infrared (IR) divergent 
one-loop contributions to the n-parton configuration, which are the virtual corrections, and 
the NLO tree level matrix elements with (n + 1) partons, the real corrections. The tree- 
level matrix elements have to be integrated over the phase space of the additional parton, 
which gives rise to collinear and soft singularities. After renormalization, the singularities 
in the virtual and soft / collinear contributions cancel and remaining poles are absorbed into 
parton distribution functions. One wants to integrate most of the phase space of the real 
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corrections numerically, but one needs to find a procedure to calculate the soft /collinear 
contributions analytically as to explicitly cancel the poles from the virtual corrections. The 
two basic methods to perform these integrations are the subtraction method |I5|-P7] and the 
phase-space slicing (PSS) method lfL8 -|2li (see also |22] for a review). 



In the following we will make use of the PSS method and therefore discuss this method 
further. To illustrate the method, we rely on the classical example given by Kunszt and 
Soper ]I6 |. We label the LO Born contribution as a LO = a B . The NLO cross section is 
given by the sum of the Born cross section and the virtual and real corrections, a v and a R : 

+ a v + a R = a B + C v - lim ^(0) + C —Fix) . (1) 



a NLO = a B 



Here, F(x) is the known, but complicated function representing the (n + l)-parton matrix 
elements. The variable x represents an angle between two partons or the energy of a gluon, 
the integral represents the phase-space intergation that has to be performed over the addi- 
tional parton. The singularity of the real corrections at x — > is compensated by the virtual 
corrections, given by the pole term and some constant, Cy. In the PSS method, the integral 
over the real corrections is divided into two parts, < x < 8 and 8 < x < 1. We note that 
the technical cut-off 8 should lie within the n jet region, i.e., if we define y cut = s m i n /Q 2 , 
then we should have 8 < y cut . If the cut-off parameter is sufficiently small, 8 <C y cnt < 1, 
one can write 

a R = f 1 dx F{x 
Jo X 

1 dx „, , , _ ,. 1 



/ —Fix) + F(0) ln(5) + lim -F(0) , (2) 

J s x e-^o e 



where the integral has been regularized by the term ggested by dimensional regu- 

larization. The pole is now explicit and the NLO cross section finite: 



C V + f 1 —F(x) + F(0) ln(5) . (3) 

J5 X 

Clearly, the real corrections a R should not depend on 8, and the logarithmic 8 dependence 
of the last term in eqn @ should be canceled by the integral, which sometimes is numerically 
difficult for very small parameters 8. However, an improvement of the above solution is 
possible by using a hybrid of the PSS and the subtraction methods, suggested by Glover 
and Sutton [P^j . In this method, one adds and subtracts only the universal soft /collinear 
approximations for x < 8, such that 

a R = lim { f 1 —x*F{x) - F(0) f —x* + F{0) f —A 
[y x Jo x Jo x J 

-/ —F{x)+ —F(x)-F(0) +F(0)ln(5) +lim-F(0) . (4) 

J8 X JO X _ _ e ^0 6 

A cancellation between the analytical and numerical terms still occurs, however only the 
phase space is approximated, so that this method is valid at larger values of 8. In the 
case where the phase-space is not approximated for small x the hybrid method becomes 
independent from 8. 
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B. Event generation with positive NLO weights 



We would now like to improve the above NLO jet cross section by including the PS 
in the n-jet region. We keep the jet definition, which separates the available phase space 
into two complementary regions, namely the n-jet region for Sy < s m i n and the (n + l)-jet 
(hard) region for > s m i n . The NLO corrections, i.e., the virtual and the soft /collinear 
corrections, will occur in the n-jet region, so that the n-jet exclusive cross section has a 
NLO normalization. On the other hand, the hard (n + l)-jet region contains tree-level 
contributions only. 

In current event generators, a weight will be generated in the n-jet region by calculating 
the LO, tree-leel matrix element and then redistributing this weight with help of the PS 
algorithm. The same jet criterion used to separate the matrix element and PS regions is 
used to veto events from the PS which would lie outside the n-jet region, as to avoid double 
counting. The (n + l)-jet region is described by the hard (n + l)-parton matrix elements. 

If we take the procedure for calculating the NLO corrections as described in the previous 
section, we could in principle calculate the weight in the n jet region in NLO. However, 
we would like to avoid the generation of large positive and negative weights, since this 
is not a practical starting point for an event generator, especially if one wishes to add 
hadronization after the showering. Therefore we suggest to choose the PSS parameter 6 
such that the weights are always positive. To keep most of the advantages of the NLO 
calculation, we aim to find a cut-off function 5 nlo (/x 2 ) that provides a cut-off parameter for 
any given renormalization and factorization scale which lies inside the n jet region, such that 
the sum of the Born, soft /collinear and virtual contributions are exactly zero. The NLO 
corrections inside the jet are then completely enclosed in the (n + l)-parton hard matrix 
elements, integrated down to the cut-off function. Only the tree-level matrix elements will 
serve as starting points for the PS. The reduced scale and scheme dependence of the NLO 
cross sections is subsumed into the scale dependent cut-off function. Staying in the simple 
example from eqn ([!]), the aimed cut-off function reads 



which is exactly what we are looking for. The NLO corrections are given completely by 
the hard (n + l)-parton tree- level matrix elements, which are positive definite and can be 
combined with the PS in the usual way. The function <5 nl ° will depend on the kinematics of 
the n parton configuration, as well as on the renormalization and factorization scales. Of 
course, one has to ensure in the event generation process that 5 nl ° < y CVLt , i.e. that the cut-off 
lies within the n-jet region. 

In general, the function (j5|) could yield a relatively large 5, so that the simple PSS method 




(5) 



Substituting 5 in eqn ([|) by 5 nl ° results in 




(6) 



s nlo 
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fails to give the correct answer for the NLO cross section.]] We can however calculate 
the missing parts of the approximated matrix elements, which is the cause of the error, 
numerically by use of the hybrid method, described above. Improving the result (H) with 
the hybrid method, one finds 



a NLO ~ 



dx , > f dx 

— F{x + / — 

X J X 

gnlo 



F{x) - F(0) 



(7) 



which gives additional terms of order S nl ° ln(5 nl °). Here, the difference between the approx- 
imate function F(x)\ x= o and the full expression is evaluated numerically. This is similar to 



the what Sjostrand and Friberg suggested in [|T^] as a starting point for the event generation. 

To summarize, we propose in close analogy to [[13] the following steps to produce an 
event in the PS region with positive weight and NLO normalization: 



Define the n-jet region by 



&ij ^ VcutQ 



(8) 



for all partons i, j. The (n + l)-jet region given by Sij > y cnt is described by the hard 
(n + l)-parton matrix elements. 

For each n-jet phase-space point calculate 5 nl ° according to eqn (H). The cross section 
for the n-jet region is then given by 



Vcut 



„NLO B i V I fi 

O n-jet =<? + O + O- 



<5 nl ° 



dx 



X 



rnlo 



F(x) + 



dx 
x 



F(x) - F(0) 



(9) 



We can write this symbolically as 

rnkn / F( X ) 



NLO 
n—jet 



5 1 



x 



5 



nlo / F(x) - F(oy 



x 



(10) 



as to express the Monte-Carlo nature of the procedure. 

In the n-jet region generate an Xk = s^/Q 2 with x fc 6 [0, y cut ] (the index k refers to the 
kth event). The weight which will be redistributed by the parton shower algorithm is 
given by 



W k = (y cnt - 5 



nlo\ 



Fix) 



x 



1 As we will see in the second part of the paper, for DIS single-jet production at 0(a s ) the cut-off 
5 nl ° actually is sufficiently small for the soft and collinear approximations to be valid. 
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for x k E [5 nlo ,y C ut], where W > by construction, or by 

\ X J x =x k 

for Xk G [0, <5 nl °], where W will be rather small. 

• Use the same jet resolution criterion as in the first step to veto events from the PS 
that lie within the (n + l)-jet region. 

III. INCLUSIVE SINGLE-JET PRODUCTION IN DIS 

In the following we apply the ideas of the previous section to inclusive single-jet produc- 
tion in DIS eP-scattering at 0(a s ). For this, we recall the formulae for calculating single-jet 
cross sections at NLO with the standard PSS method and describe the modified version of 
the PSS method, which is used to evaluate the cut-off function that can easily be used in 
existing event generators for electroproduction. We numerically compare the modified with 
the standard method. 



A. Single-jet cross section up to 0(a s ) 

In eP-scattering 

e(k) + P(p) -> e(fc') + X (11) 

the final state with a single jet is the most basic event with a large transverse energy E T in 
the laboratory frame. The lowest order O(a s ) partonic contribution to the single-jet cross 
section arises from the quark parton model (QPM) subprocess 

e(k) + q(po) -> e(fc') + q( Pl ) (12) 

and the corresponding anti-quark process with q <-> q. The partonic cross section for this 
process is given by 



where 



and 



\M q ^ q \ 2 = 32 



d£ q = a e 2 q \M q ^ q \ 2 (13) 



1 1 ( 4 ™) 2 flA^ 

(14) 



(p .k) 2 + (p .k'Y =8s 2 (l + (l-y) 2 ) (15) 



Here, s = xs, with s — (k + p) 2 , denotes the partonic center of mass energy squared, a is 
the electromagnetic coupling and Q 2 = 2{k.k') is the photon virtuality. The total DIS cross 
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section can be directly obtained by integrating out the complete phase space of the final 
state parton. 

At NLO, the single-jet cross section receives contributions from the real and the one- 
loop virtual corrections. The real corrections consist of the photon-gluon fusion and the 
QCD-Compton processes 

e{k) + q(po) -> e(k') + q( Pl ) + g(p 2 ) , (16) 
e{k) + g(p ) -> e(fc') + q{ Pl ) + q(p 2 ) , (17) 

together with corresponding anti-quark processes. The partonic cross sections for one-photon 
exchange are 

a s tf° qg = e 2 q (iira a (n R )) \M q ^ qg \ 2 , (18) 
otaCr^Zqq = e 2 q (Ana a {fji R )) \M g ^ qq \ 2 , (19) 

with 

3 (Pl-P2)(P0-P2) 

\M g ^ qq \ 2 = -- \M q ^ qg \ 2 (p <-> -p 2 ) 

= 16 C*-*0 (W + ^) 2 + ^) a + ^) a . (21) 

Color factors (including the initial state color average) are included in the squared matrix 
elements. Note that the initial state spin average factors are included in the definition of <7o 
in eqn (|T5|) and that the results in eqns (|20|j2T|) contain the full polarization dependence of 



the virtual boson. 

As already discussed in the previous section, the real corrections inherit characteristic 
divergencies, which are the initial and final state soft and collinear singularities. These 
can be separated from the hard phase-space regions by introducing a cut-off parameter 
s m i n ||l8l-|22|j . The hard part can be integrated numerically, whereas the soft / collinear part is 
treated analytically. The analytical integrals can be performed in n = 4— 2e dimensions. The 
poles which appear in e cancel against poles from the one-loop corrections. Remaining poles 
in the initial state are proportional to the Altarelli-Parisi splitting functions and are absorbed 
into the parton distribution functions (PDF's) of the proton, fi(x,fip) for i = q,q,g. UV 
divergencies in the one-loop corrections are absorbed into the running coupling constant 

The 0{ot s ) corrections to the O(a s ) Born term are known for quite some time |24| and 



the one-jet inclusive final states have been discussed in [[H],[^J • Since we will later on rely on 



programs provided together with the MEPJET Monte Carlo |26| for tabulating the integrals 
that occur for the initial state corrections, we will here apply the method of crossing functions 
as used in 0] and outline in [27], which is fully equivalent to the results in |^T,25|. We take 



over the notation in [57]. The finite part of the NLO partonic cross section, which is arrived 
at by summing up the virtual contributions and the singular parts of the two-parton final 
state is given by the expression 
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a, a, 



NLO 
q-^q 



a s <T e 2 \M q ^g\ 2 fC, 



Q 2 



(22) 



The finite parts of the virtual corrections factorize the Born matrix element. The factor 
JCq-,q, depending on both s min and the invariant mass of the hard partons 2p Q .pi = Q 2 , is 



given by 
K„ 



-q~*q\ 



Nc 
2tt 



-In 2 




* 2 1 ^ 



(23) 



where Nc = 3 is the number of colours. K, q ^ q may be crossed in exactly the same manner 
as the usual tree level crossing from the K, factor in e + e~ — > 2 partons as given in eqn (4.31) 
with n = in Ref f2| or in eqn (3.1.68) of p2| . Thus, eqn (|23"D includes also the crossing of 
a pair of collinear partons with an invariant mass smaller than s m i n from the final state to 
the initial state. This 'wrong' contribution is replaced by the correct collinear initial state 
configuration by adding the appropriate crossing function contribution to the hadronic cross 
section, which takes also into account the corresponding factorization of the initial state 
singularities, encoded in the crossing functions Cl^ s for valence and sea quark distributions. 
The crossing functions for an initial state parton a, which participates in the hard scattering 
process, can be written in the form 



with 



^>min) 



Nc 
2tt 



A a {x,n F ) In 



\x 2 F 



(24) 



A a (x,fx F ) = ^2A p ^ a (x,n F ) 



(25) 



and 



B^(x,^ F ) = ^B^ a (x^ F ) 



(26) 



The sum runs over p = q,q,g. The individual functions A p _> a {x, fj, F ) ana ts^ a {x, fi F ) are 
stated in the appendix. In particular all plus prescriptions associated with the factorization 
of the initial state collinear divergencies are absorbed in the crossing functions C^ 13 which 
is very useful for a Monte Carlo approach. We note that although the two-parton final state 
contributions (|20D and (pip contain the full polarization dependence of the virtual photon, 
the singular contributions occur only for the transverse photon polarization. 

Taking into account now virtual, initial and final state corrections we can write the 
hadronic cross section for the one-parton final state up to 0(a s ) as 



lparton / \ 
^had l s minj 



do £ e 2 J dx dPS ik ' +1) [/,(£, fi F ) (l + a s (fi R ) )C q ^ q (s miQ , Q 2 



i=q,q 



\M q ^q\ 2 



(27) 



To obtain the final, s m i n independent result, one also has to add the contribution containing 
the two parton final state, integrated over those phase-space regions, where any pair of 
partons i,j with Sy = (pi + Pj) 2 has > s m i n : 
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cr h^d rt ° n ( s min) = °"o E e i / dx dPS {k ' +2) 4ira s (/j, R ) 



i=q,q 



\Sij\>S 



'iK'mm 



+ \fg(x,H F ) \M, 



(28) 



The Lorentz-invariant phase space measure dPS^ k contains both the scattered electron 
and the partons from the photon-parton scattering process and is defined as 



(29) 



The bremsstahlung contribution in ( p8|) grows with In 2 s min and In s min with decreasing 
s m ; n . Once s m ; n is small enough for the soft and collinear approximations to be valid, this 
logarithmic growth is exactly canceled by the explicit — In 2 s min and — In s min terms in K, q ^ q 
and the s m j n dependence in the crossing functions. 



B. Cut-off function 



We are now in the position to reformulate the PSS method for the single-jet inclusive 
cross section for our purposes. As explained in section 2, we wish to avoid the NLO one- 
parton contributions contained in eqn ( p7|) completely. Integrating out the delta-function 
( P5p in eqn (B^J) we obtain, omitting scale dependences, 



j lparton 



27TO; 



dxdQ 2 xQ 4 



1 + (1 - yf) £ e 2 x \Mx) (l + a s JC q ^ q (Q 2 )) + a s Cf^x) 



(30) 



i=q,q 



The s m i n -dependence of this one-parton cross section is canceled by the respective unresolved 
two-parton cross section for each phase-space point (x, Q 2 ). In order to avoid the one-parton 
final states, it will be sufficient to chose an appropriate value of the cut-off parameter (which 
we denote as s^Jfn) for each phase-space point (x, Q 2 ), so that 



i lparton 
aa ha,d I nlo 



dx dQ 



2 V^mm; 



. 



To solve eqn (|31f) for s m i n , it is sufficient to solve the equation 

Y e 2 fi(x, n F ) (l + a s (fi R ) IC q ^ q (s min , Q 2 )) + a s (fj, R ) C 4 MS (x, fi F , s r 



i=q,q 



. 



(31) 



(32) 



The s min dependence of JC q ^ q can be seen in eqn (|23|), whereas the s min dependence of Cj 



is given in eqn (|24|) . For convenience, we define the sums 



F = 


Ee 2 


fi(x,l^ F ) , 




t=q,q 




A = 




Ai(x,n F ) , 




i=q,q 




B = 


Ee 2 






i=q,q 





(33) 
(34) 
(35) 
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and the functions 



V 

4 



In 



M 2 



3 9 A 

4 + 16 F 



In 2 




7T" 

y 



N c a s F F I M 2 



(36) 
(37) 



which are independent of s min up to 0(s min ). We have introduced some arbitrary scale M 2 
to keep the functions rj and ip dimensionless. The solution of eqn ( |3"2D is then given by the 
solution of the quadratic equation 



In 2 (^]-2r,ln(^) =V- 



\M 2 ) 



We find for s" 1 - 



nlo 
min 



"'min 



cxp 



ln(M 2 ) + r] - Jt] 2 + 4> 



(38) 



(39) 



where we have taken the smaller of the two solutions, since we require s m j n to be sufficiently 
small for the soft and collinear approximations to be valid. The ln(M 2 ) dependence in ([39]) 
cancels in the sum of the individual terms in the exponent. 

Inserting the s^ m function into eqn (p8| ) as a lower integration boundary for each phase 
space point (x, Q 2 ) will give the complete answer for the single-jet cross section in NLO. 
This is well suited for the purpose of combining matrix elements in NLO with the PS. It is 
important to note that the s m ' m function depends on the factorization and renormalization 
scales, so that the improved scale dependence of the NLO cross section is preserved in our 
modified approach. A crucial point, which we will study in detail in the next section, is 
whether the s^ m function obtained with eqn ( ft5|) is small enough for the soft and collinear 
approximations, made to evaluate the expressions fl23|) and (|24j), to be valid. 



C. Numerical results 

In this section we numerically investigate the solution (^). We look at the size of sj^ m 
for given x and Q 2 and study the effect of scale changes on s^f n . Furthermore, we check 
whether NLO single-jet inclusive cross sections obtained by integrating out the two-parton 
contributions down to s^|° n gives the same result as in the conventional approach, where 
one-parton and two-parton contributions, separated by some fixed s min , are summed. 

We start by looking at the sj^° n function in the region given by x G [10 -4 , 10 _1 ] and 
Q 2 e [10, 10 4 ] GeV 2 . We produce all results for one-photon exchange, i.e., neglecting possible 
contributions from Z-exchange. We employ the MRST [28] parton distributions for the 



proton and use the integration package provided with MEPJET to calculate and tabulate 



the crossing functions [26 for these parton distributions. This makes it numerically very 



convenient to use the function s^|° n , eqn fl31f).p] In Fig. |I]we have plotted s^f n as a function of 



2 The FORTRAN code for the s^f n function can be obtained upon request from the author. 
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Q 2 for the four fixed values x = 10~ 4 , 10~ 3 , 10~ 2 and 10 _1 for the scales \i = Hr = Hf = (,Q 2 
with £ = -7,1 and 4. We find values around 2 GeV 2 in the small Q 2 region, whereas they 
rise up to values between 100 and 200 GeV 2 for the largest Q 2 values. The s min values are 
larger for smaller x. The scale variation leads to small changes of the s m i n values. The 
scale variation in the actual cross sections will be still smaller, since the s m i n dependence of 
the cross sections is logarithmic. For the two larger x values, the smaller scales leads to a 
larger value of s^° n which will therefore produce smaller cross sections. For the two smaller 
x values there seems to be a compensation between the renormalization and factorization 



scale variations, leading to a very small overall variation in s^° n , especially at large Q . 

Next, we numerically compare the standard PSS method with our modified approach. 
The following comparisons are done for HERA conditions, i.e., E e = 27.5 GeV and E p = 
820 GeV, giving ^ = 300 GeV. A cut of E' e > 10 GeV is applied to the final state 
electron and we choose y £ [0.04, 1]. We take the same Q 2 region as above, namely Q 2 £ 
[10, 10 4 ] GeV 2 . Jets are defined in the laboratory frame with the kr algorithm with E^ h > 
5 GeV and |f/ ab | < 2. All cuts together restrict the x range to be x £ [10~ 3 , 1]. The 
numerical results for the one-jet inclusive cross sections in the following are produced with 
MEPJET f26|. 

In Fig. |2| we plot the NLO cross sections for the one-parton final states, which include 
the Born term, the virtual corrections and the soft and collinear contributions, together 
with the hard two-parton final states and their sum as a function of s m j n for four different 
Q 2 regions, integrated over the whole x-range. In Fig. ^| a we see how the logarithmic s m ; n 
dependence of the two-parton final states is compensated by the one-parton final states to 
give an s min independent result of a = 11.66 ± 0.02 nb up to values of s min ~ 10 GeV 2 . 
Above that value a slight variation of the sum can be observed and the s m j n independence 
is no longer ensured. For s min > 30 GeV 2 the one-parton final state obviously fails to 
give a correct s min dependence and the sum of one- and two-parton final states strongly 
decreases. We note that at even larger s min values, the one-parton final states will again 
give zero, which is the second solution of fl3~2|) which we rejected in (|39|). As an important 
result, one also sees that the value of s m i n , for which the one-parton final states vanish and 
the two-parton final states give the full answer is well within the s m i n independent region. 
Indeed, after we have introduced our s^° n -function into the MEPJET program we found that 
the one-parton final states did give zero and as a result for the two-parton final state we 
found er = 11.59 ± 0.01 nb, which agrees with the answer for small s m i n very well. Similar 
results hold for the larger Q 2 ranges, Fig. [2] b-d. The point at which the NLO one-parton 
final state contributions vanish are well within the s m j n independent region. This holds also 
for the largest Q 2 values, where the absolute size of the s m i n function is rather large of the 
order of 100 GeV 2 , as we have seen in Fig. [I]. At the largest Q 2 values the results seem to 
become even more stable with respect to the s m i n dependence. 

For a more detailed study we have calculated the single-jet inclusive cross section for 
nine different bins in x and Q 2 , namely Q 2 £ [10, 10 2 ] GeV 2 , Q 2 £ [10 2 , 10 3 ] GeV 2 and 
Q 2 £ [10 3 ,10 4 ] GeV 2 together with x £ [10" 3 ,10" 2 ], x £ [10" 2 ,0.1] and finally x £ [0.1,1]. 
The actual bins are summarized in Tab. |. In addition, we have tested the scale dependence 
by varying the squared renormalization and factorization scales together by a factor of 4, 
i.e., /i 2 = ji 2 R = fip = £Q 2 with £ = |, 1,4. The results are shown in Tab. 0-03] in pb, also 
indicating the relative difference A = |er st d — c"mod | /o"std of the standard PSS method to the 
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modified PSS. For all Q 2 intervals we find agreement of our modified approach compared 
to the s m j n independent standard approach to around one percent or better. The overall 
scale dependence is small, indicating a very good perturbative stability, as to be expected. 
However, it was not our intention to test the scale dependence, but to test whether our s^|° n 
function would reproduce the scale behaviour correctly. 

This concludes our numerical studies, showing the equivalence of the standard PSS 
method with our modified approach by integrating out only the two-parton final states 
down to a dynamical s^° n function given by eqn (|39|). 



IV. SUMMARY AND OUTLOOK 

We have given a prescription for combining fixed NLO matrix elements with PS's within 
the PSS method. It consists in removing the Born, virtual and soft /collinear contributions 
from the NLO cross section for the n-jet region by adjusting the PSS parameter s min for 
each phase space point and each scale These contributions are then included in the 

hard part of the NLO matrix elements, which are positive definite. This allows to directly 
redistribute the weight provided from these matrix elements with a PS algorithm. 

For the case of inclusive single-jet production in eP-scattering at 0(a s ) we have cal- 
culated the dynamical s m ; n parameter for each phase space point x and Q 2 and each scale 
fiR,fiF- We have numerically compared the standard calculation with our new approach of 
evaluating fixed NLO contributions and found the new approach to give reliable results. We 
especially found that the values of s m i n , for which the virtual plus soft /collinear contribu- 
tions vanish, are small enough for the soft and collinear approximations, used in in the PSS 
method, to be valid. We note that the cut-off function has been successfully implemented 



in the RAPGAP event generator p9 |, which includes the 0{a s ) tree level matrix elements. 
Numerical results and comparison to data will be discussed in a forthcoming paper. 

The next, more complicated step is the case of dijet production in eP-scattering, which 
is especially interesting because it allows a precise determination of ot s or the gluon density 
in the proton. NLO calculations in the PSS method from which the cut-off function can 



be determined are available pl| , p6| , p7| , |30| . It might turn out that the PSS method has to 
be suplemented with the hybrid method to numerically evaluate terms of order s min ln(s m i n ) 
as outlined in section 2. NLO calculations within the subtraction method are available for 
dijet production in eP-scattering [|17[], so that the expressions needed for the hybrid of PSS 
and subtraction method can readily be evaluated. We finally note that the 0(a 2 ) tree level 
matrix elements for eP-scattering interfaced with the PS are not yet available in a working 
Monte Carlo event generator. 
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APPENDIX A: CROSSING FUNCTIONS 



For reasons of completeness, in this appendix we collect from pTj the definitions of the 



functions A p ^ a and B^ a which are needed to compute the crossing functions . The 
functions are defined via a one dimensional integration over the parton densities f p , which 
also involves the integration over () + prescriptions. The finite, scheme independent functions 
A p _ a (x,fx F ) are given by: 

1 dz „ , , , f (llJVo - 2n f ) 



4™ = J x - f g (x/z,ii P ) | — 5(1 - z) 

+ 2 ( 5 ^- + ^ + ,(l-,))} (Al) 

^-/^^-^{^-^KiS)} (A2) 

4™= / -SMWf) 7 (A3) 

J x Z 4 

4™ = / -M*/z,Hf) t P&(z) (A4) 

J x Z 4 

The scheme dependent functions B¥^ h (x, fip) are given by: 



5 



A z J9K ' ^ J \\3 18 9Ac/ V (!-^) 



-2 l v - y + z(l - z) j ln(l - z) \ (A5) 



nMS 



nMS 



1-Z) 

1 dz 



f a WWF) \ {P^\{z)Hl ~ *) ~ H e i q {z)} (A7) 
[ Y j ln(l - z) - i*> 9 (*)} (A8) 



Here, rif denotes the number of flavors and Nc = 3 is the number of colors. The Altarelli- 
Parisi kernels in the previous equations are defined by: 

P^iz) = P£$\z) = 4 + l^l + z(l- *)) (A9) 

= \p!&\*) = y ( 1 + {1 ~ z)2 ~ «) ( A1 °) 

^ 4) (*) = \p^\z) = I { Z2 + { \~_f~ e ) (AH) 



P^?\z) = \P^\*) = y (^7 - <1 - *) ) (A12) 
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The Pjj are the e dimensional part of these n— dimensional splitting functions 

P^ 9 {z) = Ip£ 9 {z) = (A13) 

P^z) = \p£Uz) = ~\ z{l - z) (A14) 

The ()_|_ prescriptions in these equations are defined for an arbitrary test function G(z) 
(which is well behaved at z — 1) as 

1 1 X 



1 

X 



dzF + {z)G{z) = J dzF{z)[G{z) - G(l)} + G(l) J dzF{z) (A15) 



The structure and use of the crossing functions are completely analog to the usual parton 
distribution function. 

The numerical integrations have been performed in a computer program, which is pro- 
vided together with the fixed order Monte Carlo program MEPJET [26|| . The results for A p ^ a 
and B^ a for different values of x and fip are stored in an array in complete analogy to the 
usual parton densities, which allows a convenient and numerically quick evaluation of the 
crossing functions. 
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TABLES 









Q 2 £ [10, IO 2 ] 


Q 2 


£ [10 2 ,10 3 ] 


Q 2 £ [10 3 ,10 


4 ] 


X £ 


[io- 1 , 


1] 






bin 3 


bin 6 




x £ 


[io- 1 , 


io- 2 ] 

J 


bin 1 




bin 4 


bin 7 




x £ 


[io- 2 , 


1U J 


bin z 




bin 5 






TABLE I. 


Bins 


in x and Q 2 (given in GeV 2 ) for 


y > 0.04. Two bins are 


kinematically excluded. 








Scale ii 2 


= Q 2 








bin 






standard PSS 




modified PSS 




A 


1 






1044 ± 3 




1037 ± 2 




0.7 


2 






5734 ± 15 




5700 ± 9 




0.6 


3 






19.99 ± 0.18 




20.22 ± 0.04 




1.2 


4 






2193 ± 3 




2214 ± 4 




1.0 


5 






1882 ± 4 




1875 ± 2 




0.4 


6 






54.10 ± 0.20 




55.60 ± 0.18 




2.8 


7 






129.3 ± 0.5 




131.1 ± 0.4 




1.4 



TABLE II. Cross sections in pb for the bins 1-7 defined in Tab. Q, comparing the standard 
PSS method with the sj^° n modified PSS method for three different scales. A gives the difference 
of both methods in percent. 
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Scale n 2 = AQ 2 

bin standard PSS modified PSS A 

1 1124 ± 4 1120 ± 2 0.4 

2 6218 ± 15 6191 ± 7 0.4 

3 19.97 ± 0.18 20.15 ± 0.04 1.0 

4 2260 ± 4 2288 ± 5 1.2 

5 1987 ± 6 1981 ± 2 0.3 

6 54.21 ± 0.22 55.69 ± 0.13 2.7 

7 130.9 ± 0.8 132.1 ± 0.2 1.0 



TABLE III. Same as Tab. |TJ for fi 2 = AQ 2 . 

Scale /i 2 = \Q 2 

bin standard PSS modified PSS A 

1 984.1 ± 3 970.7 ± 2 1.4 

2 5412 ± 14 5343 ± 6 1.3 

3 20.00 ± 0.21 20.30 ± 0.05 1.5 

4 2134 ± 5 2158 ± 6 1.1 

5 1781 ± 6 1758 ± 2 1.3 

6 54.90 ± 0.29 55.90 ± 0.21 1.8 

7 129.0 ± 0.9 130.3 ± 0.3 1.0 



TABLE IV. Same as Tab. |] for /i 2 = \Q 2 . 
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FIG. 2. Inclusive single-jet cross section for £^ b > 5 GeV, |?? lab | < 2 and Q 2 > 10 GeV 2 as a 
function of s m in- 

The 

s min values at which the one-parton contributions vanish lie well within the 
s m i n independent region for all Q 2 . 
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